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ABSTRACT 

We present results of a comprehensive asteroseismic modelling of the (3 Cephei 
variable 9 Ophiuchi. We call these studies complex asteroseismology because our 
goal is to reproduce both pulsational frequencies as well as corresponding 
values of a complex, nonadiabatic parameter, /, defined by the radiative 
flux perturbation. To this end, we apply the method of simultaneous determination 
of the spherical harmonic degree, £, of excited pulsational mode and the corresponding 
nonadiabatic / parameter from combined multicolour photometry and radial velocity 
data. Using both the OP and OPAL opacity data, we find a family of seismic models 
which reproduce the radial and dipole centroid mode frequencies, as well as the / 
parameter associated with the radial mode. Adding the nonadiabatic parameter to 
seismic modelling of the B-type main sequence pulsators yields very strong constraints 
on stellar opacities. In particular, only with one source of opacities it is possible to 
agree the empirical values of / with their theoretical counterparts. Our results for 9 
Oph point substantially to preference for the OPAL data. 
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1 INTRODUCTION 

During the last two decades, (3 Cephei stars became at- 
tractive targets for asteroseismic studies. This has started 
with the pioneering papers by Dziembowski & Jerzykiewicz 
on 16 Lac and 12 Lac (Dziembowski & Jerzykiewicz, 1996 
and 1999, respectively). The analysis of multiplets in 16 Lac 
showed that rotation rate of this star has to increase inward. 
The seismic modelling of B-type pulsators have been revived 
with the paper by Aerts, Toul, Daszyiiska et al. (2003) on 
V836 Cen, in which they also found a non-rigid rotation 
and, for the first time, an evidence for overshooting from 
a convective core in a massive main sequence star. Up to 
now, the most intensively studied /? Cephei stars have been 
u Eri and 12 Lac which were subjects of multi-site photo- 
metric and spectroscopic campaigns (Handler et al. 2004, 
Aerts et al. 2004, Jerzykiewicz et al. 2005, Handler et al. 
2006). The seismic modelling of v Eri was done by Pamy- 
atnykh, Handler & Dziembowski (2004) and Ausseloos et 
al. (2004). Stronger constraints on stellar structure can be 
obtained from asteroseismology of the hybrid B-type pul- 
sators in which two types of modes are excited simultane- 
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ously, i.e., low-order pressure and gravity modes (typical for 
/3 Cephei stars) and high-order gravity modes (typical for 
Slowly Pulsating B-type stars). Such studies were done re- 
cently by Dziembowski & Pamyatnykh (2008) for u Eri and 
12 Lac. Interesting results are coming up for another hybrid 
pulsating star 7 Peg (Handler et al. 2009). 

Until now, asteroseismic studies have relied on fitting 
only pulsation frequencies. The ultimate goal has been to 
find models which reproduce observed frequencies by chang- 
ing various parameters of model and theory, provided that 
the successful mode identification was obtained beforehand. 
However oscillation frequencies are determined mainly by 
the star's interior, thus they are rather weakly sensitive to 
the subphotospheric layers. A few years ago, Daszyhska- 
Daszkiewicz, Dziembowski & Pamyatnykh (2003 and 2005a, 
hereafter DD03 and DD05) introduced a new asteroseismic 
probe, defined as the ratio of the amplitude of the bolometric 
flux variations to the radial displacement at the photosphere 
level. This is so called / parameter; it is determined in the 
pulsation driving region. 

The / parameter exhibits a strong dependence on the 
global stellar parameters, chemical composition, opacities 
and subphotospheric convection. Therefore, a comparison 
of empirical values of the / parameter with their theoreti- 
cal counterparts can yield valuable and unique information 
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on stellar structure and evolution. The theoretical values 
of / result from linear nonadiabatic computations of stel- 
lar pulsation. To explore the asteroseismic potential of the 
/ parameter, DD03 proposed the method of simultaneous 
determination of the spherical harmonic degree, i, and the 
nonadiabatic / parameter from amplitudes and phases of 
photometric and radial velocity variations. In this way we 
can also circumvent uncertainties in mode identification re- 
sulting from pulsational input. Using this method one can 
derive also an estimate of intrinsic amplitude of pulsational 
modes. 

This method has been successfully applied to S Scuti 
and (3 Cep stars. In the case of 5 Scuti variables, strin- 
gent constraints on the subphotospheric convection were ob- 
tained (DD03 and Daszynska-Daszkiewicz et al. 2005b) . For 
the S Scuti objects studied by these authors, convective en- 
ergy transport appeared to be inefficient. For (3 Cephei stars, 
DD05 have shown that the value of / is very sensitive to the 
metallicity parameter, Z, and to the adopted opacity tables. 
Moreover, in the case of B-type pulsators, an unambiguous 
solution exists only if one combines photometry and radial 
velocity data. 

In this paper we present an asteroseismic study which 
aims at parallel fitting of the pulsational frequencies and 
corresponding values of the complex / parameter, taking 
into account instability conditions. We will call such stud- 
ies complex asteroseismology. For our analysis we chose the 
/3 Cep star 9 Oph for which time series photometric and 
spectroscopic data were acquired by Handler, Shobbrook 
& Mokgwetsi (2005, hereafter HSM05) and Briquet et al. 
(2005, hereafter B05), respectively. Seismic modelling of 9 
Oph was recently done by Briquet et al. (2007) whose com- 
putations were based only on the OP tables. Here we shall 
use both OP and OPAL opacities and incorporate the / 
parameter in the seismic model survey. 

The paper is composed as follows. In Section 2 we view 
the basic characteristics of 9 Oph. Section 3 is devoted 
to mode identification of the excited modes using two ap- 
proaches. Constraints on intrinsic mode amplitude are given 
in Section 4. Results of complex asteroseismic modelling, 
adopting the OP and OPAL tables, are given in Section 5. 
A Summary ends the paper. 



2 THE PULSATING STAR 9 OPHIUCHI 

9 Oph is the /? Cep pulsator with the B2IV spectral type and 
the brightness oi V = 3.25 mag. It is also a triple system, 
with the speckle component of the B5 spectral type (McAl- 
ister et al. 1993, HSM05). In addition, B05 detected a low 
mass companion (M < IA^q) from spectroscopic observa- 
tions; the orbital period of this component is equal to 56.71 
d. The main component of the 9 Oph system {9 Oph A) , has 
been recognized as a pulsating variable of the /3 Cep type 
already in 1922 by Henroteau (Henroteau 1922) on the basis 
of radial velocity observations. The correct value of the pe- 
riod of about 0.14 d was derived by McNamara (1957). The 
projected rotational velocity of 9 Oph A is about 30 km/s 
as measured by Abt, Levato & Grosso (2002). More infor- 
mation on the 9 Oph system can be found in the catalog of 
galactic /3 Cephei stars by Stankov & Handler (2005). A few 
years ago, the study of 9 Oph A, 9 Oph henceforth, has been 
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Figure 1. The observational error box of 6 Ophiuchi in the 
HR diagram as determined from photometry by HSM05. The 
evolutionary tracks were computed with the OP-A04 opacities, 
assuming the rotational velocity of Vrot = 30 km/s, hydrogen 
abundance, X = 0.7, two values of the metallicity parameter Z, 
and no overshooting from the convective core. We show also lines 
of constant period (0.1339 d) for the radial fundamental (n = 1) 
and the first overtone (n = 2) mode. These lines will be discussed 
later on. 



resumed thanks to dedicated photometric and spectroscopic 
campaigns organized by HSM05 and B05, respectively. 

In Fig. 1 we show the position of 9 Oph in the HR 
diagram. The evolutionary tracks from ZAMS to TAMS 
were computed with the Warsaw-New Jersey evolutionary 
code adopting the OP opacities (Seaton 2005) and the so- 
lar mixture as determined by Asplund et al. (2004, 2005, 
hereafter A04). The rotational velocity at the ZAMS was 
assumed to be equal to 30 km/s. As 9 Oph has metallic- 
ity, Z, lower than 0.015, (Daszyhska 2001, Niemczura & 
Daszynska-Daszkiewicz 2005, B05), we showed also the ef- 
fect of the Z parameter on the evolutionary tracks. The 
observational values of effective temperature and luminos- 
ity for 9 Oph were taken from HSM05. Determination of 
effective temperature from spectroscopy gives much higher 
values, log Toff = 4.36 — 4.42. For the purpose of later discus- 
sion we plotted also lines of constant period with the value 
of 0.1339 d for the radial fundamental and first overtone 
mode. 

The frequency analysis of time-series multicolour pho- 
tometry of 9 Oph, performed by HSM05, revealed seven pul- 
sational frequencies in the range from 7.1160 to 7.9734 c/d. 
The dominant frequency {vi = 7.1160 c/d) corresponds to 
the earlier detected period (0.14 d) by McNamara (1957). 
Three of these seven frequncies were identified also in spec- 
troscopy by B05. In Fig. 2, we show the oscillation spectrum 
of 9 Oph. With ph we marked peaks detected in photom- 
etry and with ph+sp those three peaks detected both in 
photometry and spectroscopy. 
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Figure 2. Oscillation spectrum of 8 Ophiuchi. Seven frequen- 
cies were extracted from photometry (HSM05), marked as ph, 
and three of them were found in spectroscopic variations (B05), 
marked as sp. 



3 IDENTIFICATION OF PULSATIONAL 
MODES OF e OPH 

For the time being, mode identification in main sequence 
pulsators cannot be obtained from oscillation spectra them- 
selves. Therefore, one has to rely on other observational 
quantities which bring information on the mode geometry. 
Such observables are, for example, amplitudes and phases 
of the multicolour photometric and radial velocity varia- 
tions. Their theoretical values are computed under the as- 
sumptions of linear oscillation and plane-parallel, temporally 
static approximation for atmospheres. These assumptions 
are justified in the case of main sequence pulsators. If the 
effects of rotation on pulsation are ignored then the expres- 
sion for the complex amplitude of the light variations in the 
passband A can be written in the form (see, e.g., Daszyhska- 
Daszkiewicz et al. 2002): 

Axil) = -l.086eYr{i,0)bUDttf + D2.i + Dlt) (1) 



where 



D2,e = {2 + £){!-£), 



Di, = - 2 



GM 



dlOB 



9cff 



(2a) 



(26) 



(2c) 



e is the mode intrinsic amplitude, i is the inclination an- 
gle and £, m are the spherical harmonic degree and the az- 
imuthal order of the mode, respectively. Symbols G, M, R, uj 
have their usual meanings. The amplitudes and phases are 
given by Ax = |^a| and (px = arg{Ax), respectively. 

The term ^ describes temperature effects, D2,e - geo- 
metrical effects, and i - the influence of pressure changes. 
The terms ^ and D3 g include the perturbation of the 
limb-darkening, fe^ is the disc-averaging factor expressed as 
the integral of the limb-darkening weighted by the Legendre 
polynomial with the I degree. Derivatives of the monochro- 
matic flux, .Fa (Toff, log g), are calculated from static atmo- 
sphere models. In general, they depend also on the metallic- 
ity parameter [m/H] and microturbulence velocity ^t. Here, 



we use the Kurucz (2004) atmosphere models and Claret 
(2000) limb-darkening coefficients. 

As was already mentioned in Introduction, the / param- 
eter describes the ratio of the bolometric flux perturbation 
to the radial displacement at the level of the photosphere: 



Re{e/yr(e,^)e---*}. 



(3) 



The value of / can be obtained from linear computations of 
stellar pulsation; because of nonadiabatic character of the 
pulsation, this quantity is complex. Here we use the pulsa- 
tional code of W. Dziembowski. 

If data on radial velocity variations are available, then 
the amplitude and phases of these variations can be included 
in a process of mode identiflcation. The radial velocity is de- 
termined from observations as the first moment of well sep- 
arated line profile, M\. According to Dziembowski (1977), 
the expression for the complex amplitude of the radial ve- 
locity variations is the following 



A\ 



eY[^{i, 0)iujR (^ug + 



GM A 



(4) 



where Uf , are another disc- averaging factors introduced 
by Dziembowski (1977). 

In principle, two approaches can be used in mode iden- 
tification. The first approach makes use of amplitude ratios 
and phases differences and employs input from the pulsa- 
tion theory, i.e., the / parameter. In the linear and non- 
rotation approximation, the amplitude ratios and phases dif- 
ferences are independent of the intrinsic mode amplitude, e, 
azimuthal order, m, and the inclination angle, i, because the 
factor ey/"(i,0) cancels out. 

The second method employs amplitudes and phases 
alone in order to determine from data the empirical values 
of / and the intrinsic amplitude dependent on inclination, 
e — eY^ii, 0), together with the mode degree, I. 

Both methods require input from models of stellar at- 
mospheres. 



3.1 Determination of the I degree using 

photometry and theoretical values of / 

Observations of 9 Oph in the uvy Stromgren photometry 
were kindly provided by Gerald Handler. From these data we 
determined amplitudes and phases of the seven frequencies 
detected by HSM05 and, following HSM05, the amplitudes 
were corrected for the light contribution from the speckle 
companion: A^"'' = 1.19A„, Al°'' = \.22A^, A'y'"' = \.22,Ay. 
Although the phase differences between photometric pass- 
bands are small, they are not negligible, as shown in Fig. 3 
for the uy pair. Therefore, we included also these quantities 
to discriminate between various £ degrees. To this end we 
use, as usually, the following discriminant: 
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(5) 



where X° , A* are observational and theoretical photometric 
observables, respectively, such as Ax^/Ax^ and ipx^ — fXi- 
N is the number of passbands and if we use only amplitude 
ratios we have to skip 2 in the denominator before the sum. 
The index Ai denotes the photometric passband to which 
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Figure 4. Identification of the £ degrees of the pulsational modes detected in the photometric variations using amphtudc ratios and 
phase differences and theoretical / parameter (x|n vs. I). These results correspond to the model with the following stellar parameters: 
M = 9.5Mq, logTgff = 4.363, Z = 0.015, and computed with the OP-A4 opacity. The Kurucz atmosphere models were adopted. The 
effect of the microturbulence velocity, ^t, on the determination of £ is shown. 
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(7) 



We perform mode identification for the whole range of al- 
lowed stellar parameters and different values of metallicity. 

In Fig. 4, we show results of our mode identification for 
the seven pulsational frequencies of 9 Oph. This is an exam- 
ple for the model with parameters: M — 9.5AIq, logTcff = 
4.363, X = 0.7, Z = 0.015, computed with the OP-A04 
opacities. The character of XtW ioi the other models is 
qualitatively the same. We consider stellar atmospheres with 
two values of the microturbulence velocity: = 2 and 8 
km/s. In some cases a value of Xt 1^^^ than 1 which 
comes from large observational errors. 

Our identifications of the £ degrees agree with those 
obtained by HSM05 who used only photometric amplitudes. 
A discussion of the £ identification for each frequency will 
be given at the end of the next subsection. 



Figure 3. Photometric diagnostic diagram for the uy Stromgren 
passbands with positions of seven pulsational frequencies detected 
by Handler et al. 2005. 



the others, A^, are normalized and ai axe the observational 
errors of the X° quantities. Because we treat measurements 
in different passbands as independent variables, the formulae 
for errors (ai) are: 



A 



v4a. 
Ax, 



a'{AxJ ^ a^{Ax,) 



(6) 



3.2 Determination of the £ degree using empirical 
values of / 

Three pulsational frequencies {v\, V2, fa) were found in both 
photometric and spectroscopic observations of 9 Oph. By 
combining these two types of data, we can identify the I 
degree of these three frequencies and extract correspond- 
ing empirical values of the nonadiabatic / parameter. In 
our analysis we use the radial velocity as determined by the 
first moment of the SiIII4553 spectral line. These data were 
kindly provided by Mary line Briquet. In our analysis we in- 
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Figure 5. Identification of the £ degrees of the pulsational modes detected in both photometry and in spectroscopy {x\ vs. £). Various 
lines correspond to uncertainties in the stellar parameters (see text for more explanation). The Kurucz atmosphere models with the 
microturbulence velocity of ^ = 2 km/s were adopted. 



eluded only these spectra which coincide in time photometric 
observations. 

The method of simultaneous determination of £ and / 
has been described in details by DD03 and DD05, and we 
will not repeat this here. Each passband, A, yields l.h.s. of 
Equations (1) whereas data on the radial velocity yield l.h.s. 
of Equation (4). This set of linear equations is solved for a 
specified £ in order to determine two complex unknown 
quantities: e and ef. With three passbands and the radial 
velocity data we have a set of four complex equations. The 
idea is to derive the value of the nonadiabatic parameter 
/ and the intrinsic amplitude e from observations, which 
fit best the theoretical and observational amplitudes and 
phases for a given £ degree. The goodness of fit is calculated 
as 

i — l * 

where A'' is the number of passbands plus one (to account 
for the fact that we use the radial velocity measurements) 
and M is the number of parameters to be determined. Here 
A denotes either the photometric passband or radial veloc- 
ity. A° and are complex observational and theoretical 
amplitudes, respectively, and ax^ are the observational er- 
rors. Because we treat amplitudes and phases as indepen- 
dent variables, the formula for the errors is 

\axf ^a^{Ax)+AW{^x), (9) 

where A\ = |^a| and (p\ — arg{A\). 

The results of mode identification by the above method 
for the first three frequencies of 8 Oph are shown in Fig. 5. 
Various lines correspond to the uncertainties in the stellar 
parameters. We calculated Xe for the models from the center 
and from the four edges of the error box in Fig. 1. 

In Table 1 we give a summary of the £ identification 
from the two approaches. From both methods, the domi- 
nant frequency is undoubtedly a £ — 2 mode and the 1^3 
frequency is a radial mode. In the case of the 1/2 frequency, 
we were able to exclude the £ — 3 degree from combined 



multicolour photometry and radial velocity data. For the 
remaining frequencies, we can rely only on mode identifi- 
cation from photometry alone. The three frequency peaks, 
and U7, are equidistance and most probably form a 
close £ = 1 triplet. Other values of £ are also possible but 
because the disc averaging effects increase very rapidly with 
£, we regard these higher degrees as much less probable and 
assume, as HSM05, that this is the £ = 1 triplet. The £ — 
identification for U7 is of course excluded because 1/7 is too 
close to 1^3 which is the radial mode. The 1/5 frequency can 
be identified either as £=1 or 2 or 3, with a lower probabil- 
ity for the £ = 3 identification. For the dominant mode, B05 
identified also the azimuthal order, m — —1, on the basis of 
amplitudes and phases across the Silll line profile. 

Although there are some regularities in frequency spac- 
ing between ui, and V2, we are more cautious about be- 
lieving that they are components of the same quintuplet. 
The regularities suggest that may be the {£, m) = (2, +1) 
mode and 1^2 - the {£, m) = (2, +2) mode, as was assumed 
by Briquet et al. (2007) in their seismic analysis of 9 Oph. 
However, from the amplitude and phase diagrams computed 
by B05 from variations of the Silll line profile, the frequency 
1/3 = 7.3697 c/d looks rather like a m = mode, because 
the phase change is close to tt around the line center (see 
Fig. 3 of BOS). 

Thus, our complex seismic analysis we will be based 
on frequencies 1/3 and ve which are well identified centroid 
modes with £ = and £ = 1, respectively, and we postpone 
studies of rotational effect until more than one multiplet will 
be well identified. 



3.3 Determination of the radial orders 

In Fig. 1 we have plotted lines of a constant period (0.1339 
d) for fundamental (n = 1) and first overtone (n — 2) radial 
mode. These line correspond to vs = 7.4677 c/d. As we can 
see, within the allowed space of stellar parameters there are 
two possible identifications for the radial order: n = 1 or 
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Table 1 . Summary of mode identification using two methods for 
frequency peaks detected in 8 Ophiuchi. 
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n = 2. To discriminate between these two identifications, 
we shall use results from the method of DD03 (Section 3.2). 

Firstly, we can calculate values of Xb as a function of 
effective temperature for the case of n = 1 and n = 2. This is 
shown in Fig. 6. The vertical lines mark the allowed range of 
logTcff. Although the fundamental mode takes lower values 
of Xb, we cannot exclude the first overtone from this graph 
because both lines reach acceptable values of xl;- 

Another way to discriminate between the two radial or- 
ders for the ^ = mode is to compare theoretical and em- 
pirical values of the nonadiabatic / parameter (see DD05). 
The photometric and spectroscopic data for 9 Oph allowed 
to determine the values of / for the Vi mode with sufficient 
accuracy to use them to this aim. In Fig. 7 we show such 
comparison. Only models with stellar parameters within the 
error box of Fig. 1 were considered. As we can see, the n iden- 
tification is clear and is the radial fundamental mode. In 
addition, lines of constant values of the instability parame- 
ter, 77 = 0, are shown. Only models to the right of this line 
have unstable radial modes. In Fig. 7, we show also the sensi- 
tivity of the theoretical / parameter to metallicity prameter 
Z. A comparison of theoretical and empirical values of / 
gives also additional support to the result that metallicity 
of e Oph is lower than Z = 0.015. 

After fixing the n order of the radial mode as pi, 
we can try to associate the radial orders also for other fre- 
quencies. A survey of pulsational models showed that the 
£ = 1 triplet can be identified as pi, whereas the dominant 
{£, m) = (2, —1) mode, as a component of the g\ quintuplet. 



4 CONSTRAINTS ON INTRINSIC MODE 
AMPLITUDES 

Having determined the value of e = eY7"(i,0) from the 
method of DD03 in Section 3.2, we can get a lower limit 
for the value of the intrinsic mode amplitude, e, and, in the 
case of a radial mode, its exact value. The available obser- 
vations of 9 Oph allow to determine the values of e for u\ 
and with a sufficient accuracy. 

The dominant mode, u\, is unequivocally the 1 = 2 
mode. Moreover, from spectroscopic time series analysis B05 
identified its azimuthal order as m = —1. With these angu- 
lar numbers, we can calculate the intrinsic amplitude of the 
dominant mode as a function of the inclination angle. The 




\„ 

Figure 6. The x\ ^ ^ function of the effective temperature 
for the frequency V3 identified as the 1=0 mode. Two possible 
values of the radial order, n, were considered. Vertical dotted 
lines indicate the allowed range of log T^g . 



result is shown in Fig. 8. As one can see, inclination angles 
close to 85 — 90° are excluded because of the node line of 
the {£,m) = (2,-1) mode. The minimum value of the in- 
trinsic amplitude of the dominant mode is equal to 0.0098 
at i=45°, i.e. 0.98% of the stellar radius. In the case of the 
I'i frequency, which correspond to the radial mode, we esti- 
mated the exact value of e as 0.0014, i.e. 0.14% of the stellar 
radius. These estimates mean that the intrinsic amplitude 
of the dominant {£,m) = (2,-1) mode is at least 7 times 
larger than that of the radial mode. Here the open prob- 
lem of mode selection emerges and such determinations can 
bring us closer to a solution. 



5 COMPLEX ASTEROSEISMIC MODELLING 

We shall perform our seismic modelling in two steps. Firstly, 
we find a family of models which fit the two centroid fre- 
quencies of 9 Ophuichi, i.e., 1^3, which was identified as the 
£ = 0,pi mode, and ve, which is the centroid of the £ = l,pi 
triplet. In the next step, from this set of models we will se- 
lect only those which reproduce, within the observational 
errors, the empirical values of / of the radial fundamental 
mode. 



5.1 Fitting frequencies of the £ = mode and the 
£ = 1 centroid mode. 

To construct our seismic models, we use two frequencies: 
1/3 = 7,4677 (the £ = mode) and = 7.8742 c/d (the 
£ = 1 centroid mode), because the 1/1 frequency is a non- 
axisymmetric £ = 2 mode with m = —1, while frequencies 1^2 
and 1/5 do not have unambiguous identification. We searched 
for models fitting these two frequencies by changing various 
parameters: the mass, M, the effective temperature, T^s, 
the chemical abundance, {X, Z) , and the core overshooting 
parameter, Qqv. The effect of overshooting was computed 
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Figure 7. Comparison of empirical and theoretical values of the 
/ parameter for the radial mode v^. The upper panel refers to 
the fundamental mode hypothesis and the lower one, to the first 
overtone mode hypothesis. Theoretical values of / were calcu- 
lated with the OP-A04 opacities, metallicity Z = 0.020, 0.015 
and 0.012, and assuming stellar parameters within the error box. 
Lines of constant values of the instability parameter, rj, are la- 
beled as r; = 0. 



according to the new formulation of Dziembowski & Pamy- 
atnykh (2008). This new treatment takes into account not 
only the distance of the overshooting but also partial mixing 
in the overshoot layer. In our seismic modelling, we relied 
on the latest determination of the solar element mixture by 
A04 and used both the OP tables (Seaton 2005) and OPAL 
tables (Iglesias & Rogers 1996). The rotational splitting for 
the 1 = 1 triplet is about 0.1 c/d, giving the rotational 
velocity of Kot = 30 km/s for the stellar radius of about 
R = GRq ; this is consistent with the Vrot sin i value derived 
by Abt, Levato & Grosso (2002) mentioned in Section 2. We 
assumed this value of Kot in all computations. 
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Figure 8. The intrinsic amplitude of the dominant mode, i/i, 
(£, m) = (2, —1), as a function of the inclination angle. The mini- 
mum value of s{ui) is equal to about 0.01. The value of e for the 
radial mode, P3, is plotted as the dotted line. 



In Fig. 9 we show the results of our seismic modelling for 
the OP data on the Oov vs. Z plane. The instability borders 
are labeled with 77 = and shown as the thick solid and thick 
dotted lines for the I = Q and 1 = 1 modes, respectively. 
Only models above those lines are unstable and reproduce 
the observed frequencies 1/3 and vq. Moreover, we depicted 
also the lines of constant mass and effective temperature. We 
show results for two values of hydrogen abundance: X = 0.7 
(the left panel) and X = 0.75 (the right panel). 

From such data, one can easily find the relation for a 
family of seismic models for a fixed value of X in the form: 

Qfov = a • M + & ■ log Teff + c ■ Z + d. 

As we can see, models with lower Z require more core over- 
shooting. This is because for a higher value of Z the con- 
vective core is larger, implying a lower value of the over- 
shooting distance. The same result was obtained by Briquet 
et al. (2007) for 9 Oph and by Dziembowski & Pamyatnykh 
(2008) in their studies of hybrid pulsators, v Eri and 12 Lac. 
A similar dependence was also obtained for lower mass stars 
(M.-J. Goupil, private communication). On the other hand, 
at a given metallicity, for a higher value of X we need higher 
value of Qov to get pulsational instability. 

All these seismic models of 6 Oph prefers lower effective 
temperatures and lower masses. The lower effective temper- 
ature limit is logTcff ~ 4.34, thus most of these seismic 
models are outside the error box. We would like to mention 
that a lower value of effective temperature was determined 
from the lUE spectra by Niemczura& Daszyriska (2005) who 
derived logTes = 4.347 ± 0.016. The ultraviolet spectra of 
the 9 Oph system are dominated by the main component, 
i.e. the early B-type star, which emits most energy in UV 
and dominates this part of spectrum. Therefore, this lower 
value of log TeS can be more reliable and it is also more con- 
sistent with seismic models. The higher value of X require 
even lower values of logTefi. The effect of the opacity tables 
will be discussed in the next section. 
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Figure 9. The overshooting parameter as a function of metallicity for the seismic models of 9 Oph found from the fitting of the U3 
frequency (the £ = 0,pi mode) and the i/g frequency (the centroid of the £ = l,pi triplet). Results for two values of the hydrogen 
abundance, X, are shown: 0.70 (the left panel) and 0.75 (the right panel). 



5.2 Fitting the / parameter of the £ = 0,pi mode 



Determination of empirical values of / puts additional and 
unique constraints on seismic models. In that way, we are 
going one step further by a requirement of fitting simultane- 
ously pulsational frequencies and the corresponding values of 
/. The sensitivity of the theoretical / parameter to metallic- 
ity is very strong and was already shown in Fig. 7. How this 
quantity depends on the hydrogen abundance, overshoot- 
ing distance and adopted opacity tables is demonstrated in 
Fig. 10. As we can see, the strongest effect comes from opac- 
ity data. The amount of hydrogen and the overshooting dis- 
tance have a minor effect on the value of /. 



From a set of seismic models found in Section 5.1, we 
singled out only those which have the / parameter for the 
radial mode consistent with the observed values within the 
errors. Having also another source of opacity data, we stud- 
ied the effect of this input on the seismic models. In Fig. 11 
we plotted similar diagrams as in Fig. 9, but at a fixed hy- 
drogen abundance {X — 0.7) and considering both the OP 
tables (the left panel) and the OPAL tables (the right panel). 
Moreover, we superimposed areas resulting from fitting em- 
pirical values of / for the radial mode. As we can see, only 
models calculated with the OPAL tables can account for 
the empirical / parameters with a reasonable value of the 
overshooting distance. An interesting fact is that with both 
opacity tables it was possible to reproduce the real parts 
of / for lower values of Qov, whereas only the OPAL data 
gave consistency also in the imaginary parts of /. This is a 
very important and revealing result. It shows how complex 
asteroseismic studies can test data on microscopic physics. 
Another support for the OPAL data is that with these opac- 
ities we get also a better agreement between stellar param- 
eters of seismic models and those derived from photometric 
or spectroscopic observations. 
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Figure 10. A comparison of empirical and theoretical values of 
/ for the radial mode of Oph. Theoretical / parameters were 
computed for different values of the hydrogen abundance, X, the 
overshooting parameter, Oov , and for two sources of opacity data. 
The A04 solar chemical mixture was adopted and metallicity fixed 
at Z=0.015. 



6 SUMMARY 

Our analysis showed a great potential of probing param- 
eters of stellar models and microphysics by means of the 
nonadiabatic parameter /. Incorporating this quantity in 
seismic modelling provides much stronger constraints and a 
new kind of information. The adequate seismic model should 
account not only for oscillation frequencies but also for the 
empirical values of /, visibility and instability conditions. 
This kind of seismic modelling we called complex asteroseis- 
mology. The aim of the complex asteroseismic analysis is 
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Figure 11. The overshooting parameter as a function of metallicity for the seismic models of d Oph found from the fitting of the 
frequency (the i = 0,p\ mode) and the frequency (the centroid of the I = l,p\ mode). Results were obtained with the OP (the left 
panel) and the OPAL table (the right panel). The hatched areas mark the models which fit the / parameter for the radial fundamental 
mode. 



to find stellar models which fit both oscillation frequencies 
and corresponding / parameters. The value of / is deter- 
mined in subphotospheric layers which only weakly affect 
the pulsational frequencies. Therefore, these two asteroseis- 
mic tools (v, f) are complementary to each other and treat- 
ing them simultaneously can improve seismic modelling of 
any kind of pulsating stars. As was already discussed by 
DD05, for B-type pulsators, the / seismic tool probes in 
particular stellar metallicity and opacities. In this paper, we 
have demonstrated that even with data on only two well 
identified centroid modes and the empirical / parameter for 
the radial mode we can derive very strong constraints on 
stellar parameters, core overshooting and, in particular, on 
opacities. 
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Our analysis of 8 Ophiuchi showed that only with the 
OPAL table it was possible to reproduce simultaneously the 
two pulsational frequencies and the empirical value of / cor- 
responding to the radial mode. Seismic models obtained with 
the OP data require enormous amounts of core overshoot- 
ing, at least Qov ~ 0.55, which we consider not physical for a 
mass and rotational velocity appropriate for 6 Oph. Compu- 
tations with the OPAL tables allowed to achieve agreement 
between observed and theoretical values of / for a wide space 
of parameters, (M, Tefi, Oov, Z), fixed by the mode frequency 
fitting. Such complex asteroseismic studies provide a critical 
and unique test for stellar opacities and the atomic physics. 

Another important results of this paper are valuable 
constraints on the intrinsic amplitude for the dominant 
{£ — 2,m — ~1) mode and for the 1 = mode. We have 
found out that the dominant mode of 6 Oph has to be ex- 
cited with the intrinsic amplitude at least seven times larger 
than the radial mode. Although, in the framework of linear 
theory, we cannot compare these values with theoretical pre- 
dictions, such estimates can give us some guidelines to the 
mode selection mechanism in the /3 Cephei pulsators. 
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